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Figure 1: The 5-point tensor coefficient £3333 in the region of a vanishing sub-Gram determinant. Blue 
curve: conventional Passarino-Veltman reduction [18], red curve: PJFry [14, 16, 17] . Figure copyright: V. 
Yundin, 2011. 



1. Introduction 

In recent years, we worked out a tensor reduction formalism for one-loop Feynman integrals 
with more than four legs. Tensor integrals are 

with denominators cj = {k — qjY — mj + is and chords qj. 

An efficient and stable numerical implementation is of high relevance for the description of 
multi-leg final states at the LHC. The state of the art has been regularly summarized at the "Les 
Houches Workshops on Physics at TeV Colliders", see e.g. [1] or the forthcoming proceedings of 
the 2012 edition. For a broader context, see also [2-4] and references therein. 

The basics of our algebraic approach were formulated in [5-7]. In a series of papers, we 
performed the explicit reduction of 5- and 6-point tensors [8- 15]. Arbitrary internal masses and 
external virtualities are allowed. Numerical singularities arising from inverse small 5-point Gram 
determinants are avoided and those from inverse small 4-point Gram determinants are safely eval- 
uated by series of 3-point functions in higher dimensions with improvements by Pade approxi- 
mants [14]. Based on that, the numerical C-i~i- package PJFry was made available open-source 
[14, 16, 17]. An instructive example of typical output for a 6-point function is figure 1; for details 
on the kinematics see [17]. 

Since then, the interesting case of contractions of tensor integrals with external momenta was 
studied. After such contractions, the resulting scalar quantities are compact linear combinations 
of the basic scalar integrals with factorizing, simple combinations of signed minors and scalar 
products of external momenta [19, 20]. 



2 



Tensor reduction ofFeynman integrals 



Tord Riemann 



These developments were summarized at the conference, but because they were described 
already in earlier write-ups, we restrict ourselves here to the above quotations. Quite recently, we 
also studied the question of how to extend the tensor reductions to more than n = 6 external legs. 
We again relied on the earlier treatment given in [7]. The reductions may be performed consistently 
in terms of 1- to 4-point scalar functions, as was done for n <6, but one has in this framework as 
an additional element of the basis the 5-point function /j ^ ind = 6 — 2s dimensions. 

This is described in section 2. The additional integral /j''^' is finite, but it is not contained in 
packages with scalar integrals like LoopTools/FF [21, 22], QCDLoop/FF [23, 22], or OneLOop [24]. 
In section 3, we indicate how to evaluate this integral directly. 

After the conference we studied the approach due to [25, 26]. It allows to avoid the use of 
j[d+] rpj^- together with the effects of contractions with external momenta, is studied for 7- and 
8-point tensor integrals in [27] , but the details for higher tensor ranks have to be worked out yet. 



2. Tensor reduction for 7-point functions with recurrence relations 

Let us shortly repeat the reasoning for 6-point functions, and then switch to n > 6. In [12] we 
have set up the tensor reduction for the 5-point functions of rank R, expressing any (5,7?) pentagon 
by a (5,/? — 1) pentagon plus (4,7? — 1) boxes: 



'5 ~ ^5 ^0 L^^4 

.v=l 



Q'o -Lit Q^- 



(2.1) 



In [14] we have also given the corresponding homogeneous formula for 6-point functions 

6 



with p(, = and 



j^l...^R-l^l 

^6 



s=l 



(2.2) 



(2.3) 



!=1 



\QJ6 



Now we follow [7]. For a reduction of 5-point functions two recursion relations are in principle 
sufficient: 



{d-Y,Vi+m 



(d+2) 



1 

a 
1 



n k=\ 



n k=\ 



I'' 



(2.4) 
(2.5) 



where i*,j^,k* act by shifting the indices Vi,Vj,Vk by ±1. They do not work out for 6-point 
functions since ()^ = for n > 6. A further recursion relation, which does not decrease dimension, 
reads 



{d) 



[o)n k=i 



Ok 



(2.6) 



'See also [28]. 
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If reduction of dimension is acceptable, then the above double sum is often reduced by means of 

I v,j+/r^) = (2.7) 
Further one uses (4.4) of [10] together with (4.5), 

I^fO) =0. (2.8) 

!=1 V'/6 

Collecting the above, 6-point functions may be reduced. 

For 7-point functions, things are a bit more involed because we have ()7 = and Cj)^=0. 
Nevertheless one can apply recursion relations based on (2.6) and (2.7), see [7, 10]: 



o) = [» + 1 - (</ + 2.)] Q ,!f-r + t («) /S'""'. (23) 

^y^- = l» + 2-(. + 2.-)l(») /^^^^^^^^^ (2^10, 



where 



[/]' = /.^l""', (2^13) 



and the [ij]lg^ is [/j]*, but without repetition of equal indices ij. 

In [7] an idea of how to proceed for 7-point functions was formulated - details, however, were 
not given. For the 7-point vector one obtains from [5] 

(=1 

The (2.1 1) with Q^ = and {^^j = yields for / = j = A: and ;c = 2 

(:;:)/r^l.®e'-o. 

Since for the 6-point function ()g = = 0, we have from eq. (55) of [7] 

7 (Rr\ (Rr\ 

^S^''= I 7#^Sr''" + 7#^r'^ /? = anyvalueof0,...,7. (2.17) 



4 



Tensor reduction ofFeynman integrals 



Tord Riemann 



Applying in standard manner the recursion for the 5 -point function 



(2.18) 



we replace /j by integrals of the type I^j and /| ■, the dimension of which must not be reduced. 
Therefore we have to apply recursion (2.6) in the form 



(2.19) 

i.e. starting at « = 2 and increasing n step by step, we obtain the desired integrals /j and For 
« = 2 we have 

:)/--'-">C)-(:0/-(:0/- 



with 



(2.20) 



ill = ('"?)' (2.21) 



Il2 = ^h{ml). (2.22) 



2 



The problematic case is the integral /g"'^' '^ for which we can write similarly to (2. 17) 



isr)l jld+],rs (;>) 



1 (Rr\ (Rr\ 

i\d+],r 



6,! H iRr\ ^5J ~^ fRr\ ^6 

(2.23) 



The /g cannot so easily be eliminated as in the case of the 6-point function, where the vanishing 
of (2.8) was used. Inserting (2.23) into (2.17), there occurs (f/)^, i.e. quadratic, so that the right 
hand side of (2.8) does not vanish. 



3. Numerical evaluation of higher-dimensional scalar integrals 

Since in the described approach the higher dimensional integral l^^^''' cannot be eliminated, it 
is needed to investigate the possibilty for its numerical evaluation - if one wants to continue with 
this approach. First of all one would reduce the 6-point function to 5 -point functions by 

7 /Rr\ 

^r''' = l7#4'"^'^ (3-1) 
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such that the problem is shifted to the numerical evaluation of the /j . This, however, is a well 
known pathological case, since reducing the 6-dimensional 5-point function to 4-dimensional 5- 
and 4-point functions, one has 

/ON - 
\sj 5 js 



O5 



0. 



1 



d-4 



(3.2) 



i.e. for <i = 4 one meets a division by zero, ^. In [29], however, an interesting approach has been 
proposed, which can also be applied to handle this case. In fact, the idea is to go to even higher 
dimensions, which provides good numerical stability.^ One-loop «-point integrals in arbitrary di- 
mension d can be expressed in standard manner in terms of Feynman parameters as 



2 Jo Jo 



■ , ,d/2-n 

"-fn— 1 •'n f^n ; 



where 



(3.3) 



(3.4) 



and hn is a polynomial in the integration variables as well, containing also masses and momenta 
squared. The idea is to transform these integrals into integrals of higher dimension D = d — 2e + 
In — 2. For small £ the expansion of ij'f^^'^ ^' in e then reads [29] 

/i'+'"-'^=r(i + £) 

where s,, can be written as 

1 
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-Rn + 0{e) 



(3.5) 



(3.6) 



with Y: i 



- [qi — qjY + mj + mj and it is 

Rn= dx\ . . . / dXn-\ Jn hn ln(/l,j). 

Jo Jo 



(3.7) 



Obviously in such an integral no infinities occur anymore in the integrand and numerical integration 
is straightforward. It might even be possible and useful to evaluate it analytically. In our case of a 
5-point function we have n = 5, i.e. this formula applies for the integral /|'^ 
relations, we can express (3.2) as 

\2 5 



. Applying recursion 



jld+] 



{d - 2)d{d + 2) + {d- 2)d^ £ 



\0J5 



s=l 



' 



ld+],s 



(3.8) 



Here we have now the desired representation of I^^^' in terms of higher-dimensional 4-point func- 
tions, for which standard recursions can be used to reduce them to integrals in generic dimension 
d. In this way 7-point functions could finally be dealt with. 



Other approaches to the numerical treatment of higher-dimensional scalar functions are e.g. [30, 31]. 
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